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We use a variational procedure to study finite density QCD in an approximation in 
which the interaction between quarks is modelled by that induced by instantons. We find 
that uniform states with conventional chiral symmetry breaking have negative pressure with 
respect to empty space at all but the lowest densities, and are therefore unstable. This is a 
precisely defined phenomenon which motivates the basic picture of hadrons assumed in the 
£^ ■ MIT bag model, with nucleons as droplets of chiral symmetry restored phase. This suggests 

that the phase transition to a chirally symmetric phase occurs by percolation of preexisting 
droplets of the symmetric phase, and in these proceedings we expand upon our previous 
presentation of this observation. At all densities high enough that the chirally symmetric 
phase fills space, color symmetry is broken by the formation of a (qq) condensate of quark 
Cooper pairs. A plausible ordering scheme leads to a substantial gap in a Lorentz scalar 
channel involving quarks of two colors. 

m 
en 

| §1. Introduction 

00 ' 

In my lecture at YKIS97, I presented recent work done with Mark Alford and 
Frank WilczekcP on QCD at finite baryon density. This contribution to a large 
extent follows that paper, although in Section 3 I have given a significantly expanded 



Oh. treatment of chiral symmetry restoration. 

Asymptotic freedom leads us to expect that at high baryon density quarks 



behave nearly freely and form large Fermi surfaces, with interactions between the 
quasiparticles at the Fermi surfaces which become weak at asymptotically high 
density. Since the quark-quark interaction is attractive in the color 3 channel, BCS 
pairing of quarks will occur no matter how weak the interaction. Cooper pairs of 
quarks cannot be color singlets, and so a (qq) condensate inevitably breaks color 
symmetry. This breaking is analogous to the breaking of electromagnetic gauge 
invariance in superconductivity. In this color superconducting phase, the Higgs 
mechanism operates and five of the eight gluons become massive, and the U(l) 
gauge boson which remains massless is a linear combination of the photon and a 
gluon. The proposed symmetry breaking in diquark channels is of course quite 
different from chiral symmetry breaking in QCD at zero density, which occurs in 
color singlet quark-antiquark channels. 

Our goal is to explore such new forms of ordering in the context of a model 
which is definite, qualitatively reasonable, and yet sufficiently tractable that likely 
patterns of symmetry breaking and rough magnitudes of their effects can be identified 
concretely. In the course of looking for new patterns we will need to discuss the fate 
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of the old one, and here a surprise emerges: we find that a uniform phase with broken 
chiral symmetry is unstable at any nonzero density, and therefore unphysical. At 
all but the lowest densities, this instability is signalled by negative pressure, which 
presumably triggers the break-up of the uniform state into regions of high density 
separated by empty space. This suggests that in this model, baryons can be identified 
with high baryon density droplets containing massless quarks, and suggests that 
chiral symmetry restoration at high density and zero temperature can be treated 
as a percolation (or better, as we shall see, "unpercolation" ) process. In these 
proceedings, I give an alternate description of these phenomena, which has virtues 
both as a way to better understand the zero temperature physics of concern to us 
here and because it can be extended to nonzero temperatures. u We construct the 
thermodynamic potential for the model at hand, and discover that it has a first order 
transition between a state with zero baryon density and broken chiral symmetry 
(namely the vacuum) and a state with chiral symmetry restored and with a baryon 
density no greater than that of nuclear matter. Ordinary nuclear matter can then be 
identified within the model as being in the mixed phase of this first order transition: 
it consists of droplets of the high density phase separated by regions of vacuum. 
Chiral symmetry restoration occurs as the first order transition to the high density 
phase is completed as the density is increased and the droplets merge. The high 
density phase is a color superconductor, as we discuss in Section 4. t$ 

§2. Model 

We already alluded to the fact that asymptotic freedom suggests qualitatively 
new types of order at very high density. This regime has been studied If by 
approximating the interquark interactions by one gluon exchange, which is in fact 
attractive in the color antitriplet channel. Perturbative treatments cannot, by 
their nature, do full justice to a problem whose main physical interest arises at 
moderate densities. To get more insight into the phenomena, and in particular to 
make quantitative estimates, it seems appropriate to analyze a tractable, physically 
motivated model. 

We have done a variational treatment of a two-parameter class of models having 
two flavors and three colors of massless quarks. The kinetic part of the Hamiltonian 
is that for free quarks, while the interaction Hamiltonian is a slight idealization of 
the instanton vertex!^ from QCD, explicitly: 

H I = -K J d 3 x i> Rla ^ L k< i> mp iPu S e kl (35 a y s - 8 a s 5 f3 1 ) + h.c. , (2-1) 

where 1,2, k, I are flavor indices, ot,/3, 7, S are color indices, repeated indices are 
summed, and the spinor indices are not shown. The overall sign is chosen negative for 
later convenience, so that K > results in chiral symmetry breaking. Hj is not yet a 
good representation of the instanton interaction in QCD: in order to mimic the effects 
of asymptotic freedom, we must modify it in such a way that the interaction decreases 
with increasing momentum. We write Hj as a mode expansion in momentum space 
involving creation and annihilation operators and spinors, and multiply the result by 
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a product of form factors each of the form F(p) = [A 2 /(p 2 + A 2 )] U , one for each of the 
momenta of the four fermions. This factorized form is taken for later convenience, 
and is an idealization. A, of course, is some effective QCD cutoff scale, which one 
might anticipate should be in the range 300 - 1000 MeV. v parametrizes the shape 
of the form factor; we consider v = 1/2 and v = 1. Since the interaction we have 
chosen is not necessarily an accurate rendering of QCD, we will have faith only in 
conclusions that are robust with respect to the parameter choices. 

The color, flavor, and Lorentz structure of our interaction has been taken over 
directly from the instanton vertex for two-flavor QCD. This interaction properly 
reflects the chiral symmetry of QCD: axial baryon number symmetry is broken, while 
chiral SU{2) x SU(2) is respected. Color is realized as a global symmetry. There 
are other four-fermion interactions in addition to Hi which respect the unbroken 
symmetries of QCD; using Hj alone is the simplest way of breaking all symmetries 
broken by QCD, and is therefore a good starting pointrll This model of fermions 
interacting via a four-fermion interaction is one in the long line of such models 
inspired by the work of BCS as adapted to particle physics, starting with Vaks 
and Larking and Nambu and Jona-Lasinio& and studied subsequently by many 
othersEP . There is also a tradition of modelling the low-energy dynamics of QCD, and 
specifically the dynamics of chiral symmetry_breaking, with instanton interactions 
among quarks derived semi-microscopically. H This approach cannot explain the 
other main qualitative aspect of low-energy QCD dynamics, that is strict confinement 
of quarks, but it is adequate for many purposes. 

§3. Chiral Symmetry Breaking and Restoration 

We first consider symmetry breaking in the familiar pattern known for QCD at 
zero density. Working first at zero density, we choose a variational wave function of 
the form 



II cos(^(p)) + e^(P) sin(^(p))aL>)&L Q (-P) 
PA° ) 7 (3-1) 

cos(^(p)) + e^(P) edn(0*(p))ak fa (p)&t fa (-p) |0> 



with variational parameters 8 and £ depending on the modes. a\ a and b\ a create 
particles and antiparticles respectively, with flavor i and color a. This standard 
pairing form preserves the normalization of the wave function. The pairing occurs 
between particles and antiparticles with the same flavor and color but opposite 
helicity and opposite 3-momentum. It is in a Lorentz scalar, isospin singlet 
component of the chiral (2,2) representation, i.e. has the cr-field quantum numbers 
standard in this context. Following a well-trodden path we find that the energy is 
minimized when 

tan(2fl L (p)) = tan(20*(p)) = EM^K . ^ (p) + ^ (p) = n . (3 . 2 ) 

P 

For example, one could add a four-fermion interaction based on one gluon exchange. B 
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The full specification of £ R and £ L depends on spinor conventions. The gap parameter 
A is momentum independent and is defined by 

A x = 16K | o °° ^F(p) S m(9(p)) cos(0(p)) . (3-3) 
(|3-2|) and fl3-3j) are consistent only if A = or if A satisfies the gap equation 



l = m rv^P EM , (3-4) 
Jo 2vr 2 /F*(p)Al+p 2 



Note that, unlike for standard pairing at a Fermi surface, this equation does not 
have solutions for arbitrarily weak coupling, but only for couplings above a certain 
threshold value. Also note that although it is common to refer to A x as the gap 
parameter, it is best thought of as inducing an effective quark mass, which takes 
the form A x F 2 (p). In the interests of simplicity of presentation, we quote all results 
in this paper for A x = 400 MeV; we have verified that the picture we present is 
qualitatively unchanged for 300 MeV < A x < 500 MeV. Fixing A x at zero density 
fixes the magnitude (and sign) of the coupling K for each A and v. 

We now generalize this calculation to nonzero quark number density nxl It is 
most favorable energetically to fill the particle states up to some Fermi momentum 
Pf, while leaving the corresponding antiparticle states empty. That is, we replace 
|0) on the right hand side of ([3d]) by \pf), in which all particle states with [p| < pp 
are occupied. The quark number density, 

2 n . . 

n = -oPf , (3-5) 

is determined by pf but as we shall see pf is not equal to the chemical potential \x. 
The creation operators in ( |3-1| ) for modes below the Fermi surface annihilate \pf), 
and so for these modes effectively 0(p) = 0. On the other hand for the unoccupied 
states, with |p| > pp, the variational scheme is unmodified. Note that the condensate 
does not affect n, since it pairs quarks and anti-quarks. Thus we arrive at a very 
simple modification of the gap equation: the lower limit of the integral in (|3-4j ) is pf, 
rather than 0. Since the interaction is ever more effectively quenched as pf increases, 



the gap parameter A x {pp) arising as the solution of (3-4) will shrink monotonically 



and eventually vanish as pf increases. Let us define n c to be the critical density at 



which A x vanishes. 

Having obtained a definite wave function, we can evaluate the energy density in 
terms of the gap parameter. Relative to the energy of the naive vacuum state |0), 
the energy density e(n) is given by 

1 /x Pf , f°°P 2 dpp( p \ A 2 



24 v ' 16vr 2 J PF 2vr 2 2 \^ jA 2 F 4 (p) +p 2 ) 32K 



*' The reader should be warned that the next page or so is in a sense misleading; as we will see 
in due course the state we are about to construct is unphysical. 
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Fig. 1. Chiral gap (in GeV) and the pressure (in GeV 4 ) as a function of pf = (nn 2 /2) 1//3 in GeV. 
The pressure is positive at very small pf; it becomes negative at a pf which is less than 0.001. At 
n — n c , the gap vanishes and the pressure, which is still negative, has a cusp. At n — no > n c , 
the pressure crosses zero and becomes positive. 



where pjr(n) is to be obtained from ( |3-5| ). The chemical potential fi, the minimum 
energy required to add one more quark to the state, is given by \i = de/dn at fixed 
volume. We have verified that fi 2 = p 2 F + A 2 F(pp). Note that e(0) < 0, reflecting 
the fact that the physical vacuum state with its chiral condensate has a lower energy 
than the state |0). Measured relative to the physical vacuum, the energy density at 
nonzero n is e(n) = e(n) — e(0). 

The equation for A x (pp) can be solved numerically, and the resulting energy 
evaluated. Great physical interest attaches to the pressure 



P(n) 




n 



d_ 

dn 



(3-7) 



of a region with density n, where the pressure of the physical vacuum is by definition 
zero. For all values of the parameters that we consider reasonable, we find that after 
a tiny interval of very low densities at which the pressure is positive, the pressure 
becomes negative, and continues to decrease until the critical density n c . At that 
point we switch over to an essentially free quark phase, in which the energy density 
relative to the physical vacuum is 3p F /(27t 2 ) — e(0), and the pressure is given by 



P f /{2tt 2 ) + e(0). At n 



where these phases join, the pressure function is 



continuous, with a finite negative value, but it has a cusp. As density increases in 
the free-quark phase, the pressure increases monotonically, and passes through zero 
again at some density uq > n c . At n = uq, the energy per quark e/n is minimized 
and this phase is stable at zero pressure. A and P are shown for one set of parameters 
in Figure 1. 

Evidently, at all but the lowest densities (which we discuss below) in the presence 
of a chiral condensate the negative pressure associated with increasing vacuum 
energy overcompensates the increasing Fermi pressure. This negative pressure signals 
mechanical instability of the uniform chiral symmetry broken phase. There is an 
attractive physical interpretation of this phenomenon. The uniform nonzero density 
phase will break up into stable droplets of high density n = hq in which the pressure 
is zero and chiral symmetry is restored, surrounded by empty space with n = P = 0. 
Although our simple calculations do not allow us to follow the evolution and eventual 
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stabilization of the original quark cloud, it is hard to avoid identifying the droplets 
of chiral symmetric phase into which it condenses with physical nucleons. Nothing 
within the model tells us that the stable droplets have quark number 3; nucleons are 
simply the only candidates in nature which can be identified with droplets within 
which the quark density is nonzero and the chiral condensate is zero. If correct, 
this identification is very reminiscent of the MIT bag philosophy, here arising in the 
description of a sharply defined physical phenomenon. It seems quite different, 
at least superficially, from the Skyrme model, where the chiral symmetry order 
parameter changes in direction but not in magnitude within the nucleon. 

An alternative way of reformulating the results just presented is to consider the 
thermodynamic potential f2(n,A x ). Values of the gap parameter A x which solve 
the gap equation are extrema, satisfying df2([i, A x )/dA x = 0, and at these extrema 
fi/V = —P = e — /in. Using a mean field approach to evaluate 17, the following 
picture emerges. & At u = 0, 17(0, A x ) has a minimum at A x = 400 MeV which has 
P = 17 = and n = dP/d^i = 0. This is of course the vacuum solution to the gap 
equation. As /i is increased sufficiently, a second minimum of the potential develops 
at A x = 0. However, the region around the minimum describing the vacuum state 
remains unchanged, and this minimum continues to be the lowest possible 17. As 
u is increased further, the value of Q at the local minimum becomes lower and 
lower. The lowest 17 state of the system remains unchanged (eg the system is in 
the vacuum state with density zero and A x unchanged) until the chemical potential 
is increased to hq = PF( n o)- At = no, the local minimum of 17 at A x = has 
17 = 0. That is, l7(/i , 0) = = 17(^ , 400 MeV). There is now a P = 17 = chirally 
symmetric phase with density uq. Before \x can be increased any further, the system 
must undergo a first order phase transition during which the density increases from 
zero to n,Q. At all intermediate densities, the system is in a mixed phase in which 
there are regions with density no and regions with density zero. Once the transition 
is completed, u can rise further, and the lowest 17 configuration becomes one with 
17 < 0, that is with positive pressure, in which n> uq. 

What, then, of the negative pressure phase constructed previously? This solution 
to the gap equation is a local maximum of 17, with 17 > and P < 0. Were one 
to somehow construct it, it would immediately fall apart into a mixture of the two 
stable phases. The description of the preceding paragraph allows one to arrive at 
the correct description of the system (density zero, or density no, or mixture of the 
two) without ever constructing an unphysical uniform state with n and A x nonzero. 

Finally, what of the positive pressure phase at very low density? This dilute gas 
of quarks with mass A(0) has u > A(0) > /j,q. At this large value of /i, 17 has two 
nondegenerate minima. The low density gas of constituent quarks is one; the other 
is a phase with n > no, and it is this dense phase which has the lower 17. The low 
density gas of constituent quarks could only be reached if, as fi is increased above 
/xo, the phase transition to the favored high density chirally symmetric phase does 

*' Considerations similar to those we describe also lead Buballa0 to conclude that in a Nambu 
Jona-Lasinio model with an interaction which differs from the one we use, matter with broken 
chiral symmetry is unstable and nucleons can therefore only be viewed as bags within which chiral 
symmetry is restored. 
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not occur until /i is increased so high as to be above Z\ x (0). 

The picture we find is satisfying since the model finds a way to avoid having 
a phase with unconfined constituent quarks even though we have not included 
confinement in any explicit sense. This nice result is not obtained for all parameter 
values, however. For example, for v = 1 and A > 2.2 GeV, we find that (j,q > A x . 
That is, the vacuum phase becomes a gas of constituent quarks at a lower \i than 
that at which the chirally symmetric high density phase takes over. One has a mixed 
phase of droplets within which chiral symmetry is restored surrounded by regions 
not of vacuum but of regions with a nonzero density of constituent quarks. It is 
fortunate that these parameter ranges seem unreasonable. 

At a quantitative level, a naive implementation of our proposed identification 
of droplets of n = uq matter with nucleons works surprisingly well. One might 
want to identify no with the quark density at the center of baryons. Taking this to 
be two(three) times nuclear matter density yields pp(no) ~ 0.34(0.39) GeV. On the 
other hand, requiring e(no)/no to be one third the nucleon mass yields PF( n o) ~ 0.31 
GeV. Our toy model treatment cannot meet both criteria simultaneously, which is 
not surprising, but the magnitude of pf(^o) is very reasonable, ranging from 0.27 
to 0.36 for the parameters we have considered. The vacuum energy, which becomes 
the bag constant, is also reasonable in magnitude. Adding further interactions to 
Hj would obviously make a quantitative difference but the qualitative picture will 
not change, at least as long as the instanton interaction is the dominant one. 

In any case, the physical picture suggested here has significant implications for 
the phase transition, as a function of density, to restored chiral symmetry. Since the 
nucleons are regions where the symmetry is already restored, th e t ransition should 
occur by a mechanism analogous to percolation as nucleons merge. In other words, 
nuclear matter is midway through the first order phase transition described above, 
and chiral symmetry restoration is nothing more than completing this transition. 
This transition should be complete once a density characteristic of the center of 
nucleons is achieved. As n increases, the fraction of the volume of space occupied 
by the n = uq phase increases. In our toy model, the chemical potential and the 
pressure remain constant as n is increased from that of nuclear matter, which is in 
the mixed phase, to uq. This reflects the fact that we have neglected interactions 
between droplets. In nature, once these are included, some external pressure must 
be imposed in order to induce the nucleons to merge.l*j 

*' There is a related reason that interactions between droplets are needed. Although it is very 
nice to see the model favoring quarks "confined" into droplets within which they have zero mass 
rather than free as constituent quarks, it is asking too much of this toy model to expect it to 
correctly predict that nature favors droplets containing three quarks. However, augmenting the 
model to include a short range repulsion between droplets would make it unfavorable for small 
droplets to merge into bigger ones, and would mean that in the mixed phase the favored droplets 
are the smallest possible droplets. A model with color introduced only as a global symmetry cannot 
select 3-quark-droplets, but it would nevertheless be interesting to see it favoring small droplets over 
large ones. One way of introducing such an effect would be to introduce a negative surface tension 
for the interface between the n = no phase and the vacuum. With a negative surface tension, in 
the mixed phase the system will maximize the ratio of droplet surface to droplet volume, and will 
therefore consist of the smallest possible droplets of the n = no phase surrounded by vacuum. A 
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The mechanism of chiral symmetry restoration at finite density but zero 
temperature is quite different from the one we expect at finite temperature and 
zero density: it occurs by percolation among pre-formed bags of symmetric phase. 
In the toy model we have described, this transition seems perfectly smooth. As 
the density increases, one phase steadily wins over the other. In the chiral limit, 
however, this cannot be the whole story. At zero density, chiral symmetry is broken. 
At n = no when the dense phase has taken over completely, chiral symmetry is 
restored. At what density (that is at what volume ratio of the two phases in the 
mixed phase regime) is chiral symmetry restored? This question cannot be answered 
precisely using the methods we have developed to this point. It is plausible, however, 
that as soon as the density is high enough that there are no longer vacuum regions 
of infinite spatial extent, chiral symmetry is restored. Once one has only isolated 
vacuum regions, each of these isolated regions can have the chiral order parameter 
pointing in different directions in isospin space. Each isolated region will pick a 
different a direction. This means that in the system as a whole, chiral symmetry 
will be restored. In this sense, what matters is the "unpercolation" of the vacuum 
rather than the percolation of the droplets. Regardless of the precise criterion, the 
density n x at which chiral symmetry restoration occurs must be somewhat less than 
no- At densities less than n x , one has massless pions — oscillations of the orientation 
of the condensate present in the vacuum regions of infinite spatial extent associated 
with spontaneous chiral symmetry breaking. At higher densities n > n x , the pions 
need no longer be massless. It is a logical possibility that the restoration of chiral 
symmetry at n = n x is a second order transition, occuring within the "larger" first 
order transition. It is also a logical possibility that the chiral symmetry restoration 
transition is first order, between two phases with density n~ and n+, again within 
the larger first order transition between n = and n = uq. Of course, away from 
the chiral limit, with a nonzero current quark mass and a nonzero pion mass, the 
transition can be completely smooth. 



At high density, pairing of particles near the Fermi surface as in the original 
BCS scheme becomes more favorable. Our Hamiltonian supports condensation in 
quark-quark channels. The condensation is now between fermions with the same 
helicity, and the Hamiltonian selects chiral isosinglets — that is, antisymmetry in 
flavor. One can therefore have spin — antisymmetric in spin and therefore in color, 
forming a 3, or spin 1 — symmetric in spin and therefore in color, forming a 6. 

We first consider the former. A suitable trial wave function iso 



negative surface tension (or short range droplet-droplet repulsion introduced in some other way) will 
not upset the stability of either the vacuum or of the uniform phase with density n = no; it simply 
upsets the stability of a mixed phase with one large droplet relative to one with many smaller ones. 

*' Nothing misleading here; there will be no negative pressure and no droplet formation here; 
the state we are about to construct is stable. 



§4. Color Superconductivity 
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where 



G\ = II (cos(^(p)) + ^ 3 e^(P) S in(^(p))al lQ (p)aL /3 (-p)) 

fcos(0*(p)) + e ^3 e ^( P ) S m(^(p))6^ lQ (p)6^ 2/3 (-p)) (4 . 2) 

^o S (^(p))+^ 3 e^c(P) s in(^(p))a Rla (p)aT ?2/3 (-p)) 
Gjj = same, with i? <-> L 
Here, a and /3 are color indices, and we have chosen to pair quarks of the first two 



colors. 1 and 2 are flavor indices. The first term in (4-2) creates particles above 
the Fermi surface; the second creates antiparticles; the third creates holes below the 
Fermi surface. In this state, the Lorentz scalar (q ia C^q^eij £ a ^) is nonzero. This 
singles out a preferred direction in color space and breaks color SU (3) — > SU (2). The 
U(l) of electromagnetism is spontaneously broken but there is a linear combination 
of electric charge and color hypercharge under which the condensate is neutral, and 
which therefore generates an unbroken U(l) gauge symmetry. No flavor symmetries, 
not even the chiral ones, are broken. 

Note that n is now not given by ( |3-5| ) because the operators in (^^) can change 
particle number. Varying the expectation value of H — fiN in this state with respect 
to pf yields pf = M> unlike in the case of the chiral condensate. This difference 
reflects the fact that a gap in a {qq) channel does not act as an effective mass term 
in the way that a chiral gap does. Upon adding a quark, the condensate can adjust in 
such a way that the energy cost is only pp. There is, however, a true gap in the sense 
of condensed matter physics: the energy cost of making a particle-hole excitation is 
2AF 2 (/j,) at minimum. Varying with respect to all the other variational parameters 
yields £% B O + & BC = ir, Q\ B C = 6\ B C , and 

.an(20i(p» = ^ , t an(2^(p)) = ^ , ta„(2^(p)) = ^ . 

V — IX p + p n - p 

(4-3) 

Here, the gap A satisfies a self-consistency equation of the form 

] F i {p)p 2 dp f°° F 4 (p)p 2 dp 




y/F*(p)A* + {p-iJ,)* Jo ^F*(p)A' 2 + {p + nY 

r» F\p)p 2 dp 1 
Jo ^/FHp)A 2 + (p-p) 2 J ' 1 J 

The three terms in this equation arise respectively from particles above the Fermi 
surface, antiparticles, and particles below the Fermi surface. For n > the particle 
and hole integrals diverge logarithmically at the Fermi surface as A — > 0, which 
signals the possibility of condensation for arbitrarily weak attraction. 

Because the numerical coefficient in the gap equation is smaller than the 
threshold value at which one would have a nonzero A x at p = 0, A would be zero 
were it not for the logarithmically divergent contribution to the integral from the 
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Fig. 2. Gap created by the Lorentz scalar color superconductor condensate, as a function of fi = pf 
for v = 1 and (from left to right) A = 0.4, 0.8 GeV. Each curve begins where n is given by the 
appropriate no. 

region near fi. This means that at small fi, the gap must be small because the density 
of states at the Fermi surface is small. This has only formal significance, because the 
only densities of physical relevance are n = and n > uq. At intermediate densities, 
matter is in an inhomogeneous mixture of the n = and n = uq phases. As u 
increases, the density of states at the Fermi surface increases and the gap parameter 
grows. Finally, at large u the effect of the form factor F is felt, the effective coupling 
decreases, and the gap parameter goes back down. For the parameter ranges we 
have examined the gap parameter is quite substantial: ~ 50 — 150 MeV at no, and 
peaking at 100 — 200 MeV at a density somewhat higher. We plot A for two sets 
of parameters in Figure 2. The density at which the gap peaks depends on A; the 
shape of the curve depends on v\ the height of the curve is almost independent of 
both. 

As in the previous section, one can obtain expressions for the energy and the 
density, and thus derive the equation of state. We find that the equation of state is 
hardly modified from the free-quark values — the presence of the condensate typically 
lowers the energy density by several tens of MeV per fm 3 and the pressures, at equal 
baryon density, are equal to within a few per cent. This makes it very plausible that 
the color condensation makes only a small change in the density no at which a stable 
phase exists at zero pressure. We have made this argument more rigorous by doing 
a a calculation in which we consider chiral and color condensation simultaneously. 1$ 

The third color has so far been left out in the cold, but we can gain energy by 
allowing it too to condense. The available channel is the color 6. There is attraction 
in this channel, B for a pairing that is a spatial axial vector. The resulting condensate 
spontaneously breaks both color and rotational symmetry. Without modification, 
our model suggests that this axial vector gap is of order a few keV at most,EP orders 
of magnitude smaller then the scalar color breaking condensate. This is not a robust 
result, and it is worth exploring whether plausible interactions can be added to Hj 
which strengthen the axial vector condensate. 

In looking for signatures of color superconductivity in heavy ion physics and in 
neutron stars, it is unfortunate that the equation of state is almost equal to that for 
a deconfined phase with no diquark condensate. Superconducting condensates do 
modify the gauge interactions and this may have implications in heavy ion collisions. 
In our model as it stands, color is realized as a global symmetry. Breaking of this 
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symmetry generates Nambu-Goldstone bosons, formally. However, in reality color is 
of course a gauge symmetry, and the true spectrum does not contain massless scalars, 
but rather massive vectors. There is a residual SU (2) gauge symmetry, presumably 
deconfined, and there are five gluons whose mass is set by the scalar condensate. 
The scalar condensate carries electric as well as color charge. It is neutral under 
a certain combination of electrodynamic and color hypercharge. If densities above 
no are achieved at low enough temperatures that the scalar condensate forms, there 
will be a mixture of the ordinary photon and the color hypercharge gauge boson 
which is massless. Either the modification of the photon or the loss of massless 
gluons could have consequences, but dramatic effects do not yet seem apparent. 
Recently, we have described another possible signature. EJ) The formation of a color 
superconducting state lowers the energy most if the matter is up-down symmetric. 
In a heavy ion collision, in which down quarks initially exceed up quarks, as a 
color superconductor condensate forms in the densest region near the center of the 
collision it may therefore expel a few down quarks and up anti-quarks into the 
surrounding less dense regions. These will eventually become negative pions. When 
the condensate breaks up late in the collision, it likely yields equal numbers of protons 
and neutrons. The total number of protons in the final state will then be more 
than twice 79. Experimentalists should therefore look for event-by-event correlation 
between observables which suggest high density and low temperature on the one 
hand, and either an increase in the number of protons per event or an increase in 
the /tt + ratio on the other. 

Turning to neutron star phenomenology, the slowness of observed neutron star 
cooling rates may indicate that a gap in the excitation spectrum for quark matter 
might be welcome, US as this suppresses neutrino emission via weak interaction 
processes involving single thermally excited u and d quarks by exp(— Z\/T). 

§5. Discussion 

Many things were ignored in this analysis. Most important, perhaps, is the 
strange quark. Whatever the interaction(s), color superconductivity in a three flavor 
theory necessarily introduces the new feature of flavor symmetry breaking. The 
condensates considered in our two-flavor model are flavor singlets; this is impossible 
for a (qq) condensate in a three flavor theory. One particularly attractive possibility 
is condensation in the {qf C^q^e^ A e a pA) channel, with summation over A. This 
breaks flavor and color in a coordinated fashion, leaving unbroken the diagonal 
subgroup of 5C/(3) co ior x SU(3)q &vot . 

Another question concerns the postulated Hamiltonian. While there are good 
reasons to take an effective interaction of the instanton type as a starting point, 
there could well be significant corrections affecting the more delicate consequences 
such as axial vector 6 condensation. A specific, important example is to compare 
the effective interaction derived from one-gluon exchange. It turns out that this 
interaction has a similar pattern, for our purposes, to the instanton: it is very 
attractive in the a channel, attractive in the color antitriplet scalar, and neither 
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attractive nor repulsive in the color sextet axial vector. More generally, it would be 
desirable to use a renormalization group treatment to find the interactions which are 
most relevant near the Fermi surface. 

The qualitative model we have treated suggests a compelling picture both for 
the chiral restoration transition and for the color superconductivity which sets in 
at densities just beyond. It points toward future work in many directions: chiral 
symmetry restoration via the percolation (of baryons) or unpercolation (of the 
vacuum) must be further characterized; consequences in neutron star and heavy 
ion physics remain to be elucidated; the superconducting ordering patterns may 
hold further surprises, particularly as flavor becomes important. 
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